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Q\ ■ Abstract 

In this short review the role of the Hirota equation and the tau-function in the 
r*j ' theory of classical and quantum integrable systems is outlined. 
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1. Introduction. As is known, integrable models exist in two versions - classical 
and quantum. In the former one should solve equations of motion while in the latter 
the primary concern is to diagonalize operators (Hamiltonians or evolution operators). 
Usually classical integrable models admit quantizaton that preserves integrability and 
^ 1 quantum models have a well-defined classical limit, in accordance with the quantum- 

mechanical correspondence principle. 

At the same time it turns out that there are deeper links between quantum and 
classical integrable models which do not follow from the usual correspondence principle. 
Namely, it appears that classical integrable equations are built in the structure of quan- 
tum models as exact relations even at H ^ 0. And vice versa, some specific accessories 
of the quantum theory appear in solving purely classical integrable problems. There are 
several aspects of this surprising and not yet fully understood phenomenon. We discuss 
^ 1 only one - the identification of the quantum transfer matrix (the T-operator) with the 

classical r-function. 



In the main text we do not give any references. Some references with brief comments 
are collected in section 7. 

2. Quantum integrable systems: Bethe ansatz. It is noteworthy that the first 
non-trivial problem with many degrees of freedom solved exactly was quantum rather 
than classical. In 1931 H. Bethe managed to find exact wave functions of the Heisenberg 
spin chain with the Hamiltonian 
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Here cr x ,a y , <J Z are the standard Pauli matrices. The key point of the solution was a 
special ansatz for the wave function which is now called the Bethe ansatz. Later one 
or other version of this method turned out to be applicable to many other lattice or 
continuous quantum models in 1 + 1 dimensions. The typical form of the answer for the 
spectrum of the Hamiltonian is 

E = "£e(v k ). 

k 

Here e(v) is a known function (for the Heisenberg chain it is e(v) = —8 J/ (v 2 + 1)) and the 
numbers v k which are rapidities of quasiparticles are found from the system of algebraic 
equations (the Bethe equations) 



Vk + iJ i v k -vi + 2% 

In physically interesting cases, when the ground state in the thermodynamic limit is built 
by filling a false vacuum and the number of quasiparticles tends to infinity together with 
the lattice length L, the system of Bethe equations becomes an integral equation for the 
density of continuously distributed rapidities which can be solved explicitly. 

In the case of generalized spin chains with "spin variables" belonging to represen- 
tations of the group SU(N), diagonalization of Hamiltonian is performed by means of 
consecutive application of the Bethe ansatz N — 1 times (this method is usually referred 
to as nested Bethe ansatz). Correspondingly, there are several sorts of quasiparticles 
with rapidities vf\ t — 1,2, ...,7V — 1, and the Bethe equations read 

n vP-vtV+i vt^t^i vP-vt^-i _ (9) 
i v® - vt 1] - i vf - vf + 2 % vf - vl t+1) + i [ ) 

For many years a number of different versions and generalizations of the Bethe method 
were suggested. However, its main secret seems to be still unrevealed. 

3. Integrable models of classical theory: Hirota equation. It is a question about 
non-linear partial differential equations. Today we know many cases when they appear to 
be integrable. Among them are well-known examples of integrable models of classical field 
theory in 1 + 1 dimensions: sin-Gordon, the principal chiral field, sigma-models. There 
are also non-relativistic models such as the Korteweg - de Vries equation, the non-linear 
Schrodinger equation, the Kadomtsev-Petviashvili equation, the Toda chain and many 
other. They are often called soliton equations because they usually have exact solutions 
of that type (localized moving excitations that preserve their shape in the evolution). The 
integrability, i.e. the existence of infinite number of conserved quantities in involution, 
means that with each such equation one can associate an infinite hierarchy of compatible 
equations since each integral of motion generates its own evolution in time. 

The central object of the theory is a function on an infinite dimensional parameter 
space which is called r-function. In a nutshell, it provides a common solution to the 
whole hierarchy of integrable equations generated by the infinite family of Hamiltonians 
in involution. As such, it depends on infinite number of variables t = {t , t\, t 2 , t 3 , . . .} 
and satisfies infinite number of equations which can be encoded in just one functional 
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relation (the Hirota equation) 



(z 2 -z 3 ) r(t + [zi])r(t + [z 2 ] + [z 3 ]) 

+ (z 3 - Zl ) r(t + [z 2 ])r{t + [ Zl ] + [z 3 ]) (3) 

+ (zi-^ 2 ) r(t + [z 3 ])r(t + [zi] + [z 2 \) = 0. 

where Zi, z 2 , z 3 are arbitrary parameters, and we use the short-hand notation t+[z] = {io+ 
1, ti+z, t 2 + jz 2 ,t 3 + ^z 3 , . . .}. Expansion in powers of z-i yields differential equations which 
constitute the hierarchy. (To be precise, ([3J) corresponds to the modified Kadomtsev- 
Petviashvili hierarchy; in other cases the functional relations have a similar bilinear form.) 
Along with ti one may also use the variables u z which are numbered by continuous "label" 
z and which are connected with the t^s by the relations 

t = u , t k = \^£.u z z k , k>l (4) 

(let us assume that the sum is finite). In the variables u z equation becomes a difference 
equation on a 3D lattice for any triple u Zi — Ui (i — 1, 2, 3): 

(Z2-Z3) t(u 17 u 2 + 1, u 3 + 1)t(ux + 1, u 2 , u 3 ) 

+ (z 3 -z 1 )t(u 1 + 1,u 2 ,u 3 + 1)t(u 1 ,u 2 + 1,u 3 ) (5) 

+ (z 1 -z 2 )t(u 1 + 1,u 2 + 1,u 3 )t(u 1 ,u 2 ,u 3 + 1) = 0. 

The parameters z\ play the role of the lattice spacings. Similar equations can be written 
for four and more variables Ui but all of them appear to be algebraic consequences of 
equations (JSJ) written for any triple of variables. 

The sets of variables tj and u z provide complimentary descriptions of one and the 
same integrable system and transition from one to another is in some sense similar to 
the Fourier transform. Let us also note that the Hirota equation written in the forms ([3]) 
or (111) reflects a deep interrelation between continuous and discrete or difference soliton 
equations: they belong to one and the same hierarchy and turn one into another at the 
change of the independent variables (jlj) made simultaneously in the whole hierarchy. 

"In nature" r-functions or their logarithms appear as partition functions, different 
kinds of correlators and their generating functions, and effective actions as functions of 
coupling constants. 

In fact the whole variety of integrable non-linear partial differential equations can 
be encoded in one universal difference equation (J5]) for the r-function. As well as the 
equivalent equation ([3]), it is called the Hirota equation. It plays a truly vital role in 
the theory of classical (and also quantum, as we shall see soon) integrable systems. 
All known integrable equations can be obtained from it by various simple but maybe 
technically sophisticated manipulations such as continuous limit (expansion in powers of 
Zi and transition to the variables tk), imposing of reductions, choosing dependent and 
independent variables and the like. 

As any fundamental thing, the Hirota equation and closely related equation called the 
F-system often appears in different unexpected contexts. Recently it was used for finding 
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the spectrum of anomalous dimensions of composite operators in N = 4 supersymmetric 
4D Yang-Mills theory in the planar limit. 

4. Towards a synthesis of classical and quantum integrability. By quantum 
integrable system we mean, in this section, a non-homogeneous spin chain whose local 
observables ("spins" on the sites) are operators acting in finite dimensional representa- 
tions of the su(N) algebra or its g-deformation. This example is rather representative 
because many other exactly solvable quantum models can be treated, at least formally, 
as its limiting cases. 

Since the system has many or even infinitely many commuting integrals of motion 
where {J} = {Ji, J2, ■ ■ ■} is some multi-index, it is natural to diagonalize them 
simultaneously rather than one particular Hamiltonian from this family. Even better, 
one can combine them in a generating function and diagonalize this operator function. 
Schematically, it looks like 

T(t ,ti,t 2 , • • •) = E H Vh ( 6 ) 

in j 

where a^jy are properly chosen coefficients. It depends on an infinite number of auxiliary 
variables ti. Such generating function is called the master T-operator (or simply T- 
operator). It has the meaning of an evolution operator in a time determined by the 
parameters ti. The T-operator is a much more meaningful and informative object than 
any particular Hamiltonian because an important dynamical information is encoded in 
its analytic properties in the variables t{ (which are in general complex numbers). 

The commuting family of operators constructed in such a way contains, along with the 
full set of commuting Hamiltonians of the spin chain, the transfer matrices of associated 
2D lattice models of statistical mechanics (vertex models) as well as all Baxter's Q- 
operators. 

More often than not an explicit expression for the T-operator through original dyna- 
mical variables is not available but this is not an obstacle for the derivation of general 
functional relations for it. With taking into account the analytic properties, they allow 
one to solve the spectral problem for the T-operator and thus for the original Hamiltonian. 
These functional relations have a long history and are known in different forms. In a 
sense, they are at the top of the theory of quantum integrable systems. Remarkably, at 
this place the quantum theory becomes very close to the classical one. 

The key fact that provides a hint for the anticipated synthesis of classical and quan- 
tum integrability is that the most general and universal form of the functional relations 
for the T-operator T(t ,ti, . . .) (jSj) is nothing else than the classical Hirota equation 
d2J) in the variables Let us remark that since the T-operators commute for all val- 
ues of ti, the ordering ambiguity, usually accompanying quantization, does not arise 
here, and any eigenvalue of the T-operator satisfies the Hirota equation. More precisely, 
the T-operator satisfies the general bilinear relation for the r-function of the modified 
Kadomtsev-Petviashvili (mKP) hierarchy from which the equations of the Hirota type 
follow. In other words, the T-operator (or rather any of its eigenvalues) and the classical 
r-function can be identified: 

T(t , ti,t 2 , ■ ■ •) = r(t , ti, t 2 , . . .). 
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To make the statement more precise, we need some details about the T-operator. Among 
the variables t iy the first one, t , is distinguished. In order to stress this, we denote it by 
u: u — t . It is called the quantum spectral parameter. Let us consider the T-operator 
as a function of u, depending on all other t = {t\, t^, . . .} as on parameters: T = T(u, t) 
(for definiteness, one can keep in mind any of its eigenvalues). Let T\(u) be the set 
of commuting quantum transfer matrices. They are functions of u and depend on a 
Young diagram A as on discrete (multi-component) parameter. The T-operator can be 
represented as the following generating series which is the precise version of the schematic 
formula ©: 

T(M)=£s A (t)T A («). (7) 

A 

The sum is taken over all Young diagrams including the empty one and s\(t) are the 
Schur functions (which are usually regarded as symmetric functions of the variables Xi 
such that tk = r J2i X T)- The known functional relations for T\(u) imply that this object 
satisfies the classical mKP hierarchy and, in particular, the Hirota equation. The Planck 
constant h of the quantum problem is not present in the Hirota equation but appears in 
its solutions. 

In the context of the classical Hirota equation, the nested Bethe ansatz method trans- 
lates to a chain of Backlund transformations which sequentially "undress" the original 
solution to a trivial one. 

All this can be almost literally extended to quantum integrable models with super- 
symmetry. The Hirota equation for the T-operator remains the same, the only change is 
the class of relevant solutions. 

5. Connection with many-body problems of the Calogero-Moser type. An- 
other side of the classical-quantum correspondence is the link to classical integrable many- 
body problems of the Calogero-Moser type which enter into the game through dynamics 
of zeros of the T-operator. This link follows from the identification of the T-operator 
with the classical r-function and from the analytic properties of the former. 

To give a more precise statement, we need to know the analytic properties of the 
T-operator T(w,t), u = t , in the variable u (the quantum spectral parameter). This is 
one of the most important characteristics of the quantum model. In the simplest case, 
which includes finite spin chains, the T-operator should be a polynomial of u of degree 
equal to the length of the chain: 

T(u,t) = cf[(u- Uj (t)). (8) 

i=i 

Then from the identification of T(u, t) with the r-function it follows that the dynamics 
of the zeros Uj(t) in the times tj is given by the relativistic generalization of the Calogero- 
Moser model (it is usually called the Ruijsenaars-Schneider system). For example the 
equations of motion in t\ are 

Ui = J2~, — — -> iii = d tl u i . (9) 



(ui-u k ) 2 - 1 ' 



This system is integrable and has the required number of independent conserved quanti- 



ties Hj in involution. 
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Finally, let us stress the specific way of posing the classical mechanical problem for 
the Ruijsenaars- Schneider system that corresponds to diagonalization of quantum spin 
chain Hamiltonians or transfer matrices of the vertex model. The standard mechanical 
problem is: given initial coordinates and velocities of the particles Uj(0), Uj(0), find Uj(t). 
In distinction to this, in order to find eigenvalues of the T-operator one should pose the 
problem in the following non-standard way: given initial coordinates Uj = Uj(0) and 
values of all higher integrals of motion Uj, find initial velocities Uj(0). The solution 
of such a problem in general is not unique: different possible solutions correspond to 
different eigenstates of the quantum Hamiltonians. 

We also note that the Bethe equations ($Z§ can be understood as the Ruijsenaars— 
Schneider system in the discrete time t. 

6. Concluding remarks. We see that the most universal relations for classical and 
quantum integrable systems actually coincide and are given by the Hirota equation. Sum- 
ming up, we suggest the following extension of the quantum-mechanical correspondence 
principle: with any quantum integrable system one can associate a classical integrable 
dynamics in the space of its (commuting) integrals of motion. This classical dynamics 
contains complete information about spectral properties of the quantum system. 

At this point we finish the story about interrelations between classical and quantum 
integrability told "at verbal level". In our opinion, the very fact that the top points of 
the classical and quantum theories of integrable systems actually coincide takes place for 
profound reasons and requires a deeper understanding. 

7. Some references and comments. The fundamental role of the r-function in the 
theory of soliton equations was revealed in the works of the Kyoto school (see e.g. [TJ 
and references therein). The discrete Hirota equation first appeared in [2], its meaning 
for integrable hierarchies was clarified in [3] . The Bethe ansatz is reviewed in jU |5] . The 
functional relations for quantum transfer matrices in the form of determinant formulas 
were found in |HJ UJ, for the supersymmetric extension see [8]. A similarity between 
quantum transfer matrices and classical r-functions was first pointed out in [9] (see also 
|10j). where the discrete classical Hirota dynamics in the space of commuting quantum 
integrals of motion was introduced. An extension of this approach to supersymmetric 
models was suggested in [UJ . For the master T-operator see [121 02] • The dynamics of 
poles of solutions to classical integrable equations (zeros of the r-function) in connection 
with solvable many-body problems was studied in [J5J [TBI HZ]- The coincidence of the 
equations of motion for the discrete time analog of the Ruijsenaars-Schneider system and 
the Bethe equations was noticed in [18], a connection with motion of zeros of r-function 
was pointed out in |9]. For various aspects and applications of the Hirota equation (the 
T-system) and the associated F-system see review [T9] and references therein. Another 
context where r-functions of classical integrable hierarchies enter the theory of quantum 
integrable models and associated 2D lattice models of statistical mechanics is calculation 
of scalar products and partition functions with domain wall boundary conditions [20] . 
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